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We present observational constraints on a scalar-tensor gravity theory by \ 2 test for CMB 
anisotropy spectrum. We compare the WMAP temperature power spectrum with the harmonic 
attractor model, in which the scalar field has its harmonic effective potential with curvature /3 in 
the Einstein conformal frame and the theory relaxes toward Einstein gravity with time. We found 
that the present value of the scalar coupling, i.e. the present level of deviation from Einstein gravity 
(og), is bounded to be smaller than 5 x lO" 4 " 7 ' 3 (2a), and fO" 2 " 7 " (4a) for < (3 < 0.45. This 
constraint is much stronger than the bound from the solar system experiments for large f3 mod- 
els, i.e., j3 > 0.2 and 0.3 in 2a and 4a limits, respectively. Furthermore, within the framework of 
this model, the variation of the gravitational constant at the recombination epoch is constrained as 
\G{z = z rec ) - G \/Gq < 0.05(2cr), and 0.23(4<x). 
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I. INTRODUCTION 

The discovery of acoustic peaks on the anisotropy 
spectrum of the cosmic microwave background (CMB) 
opened the way to measure the early universe. The 
high quality data recently provided by the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) mission pj makes 
such a challenge feasible. The measurement of the char- 
acteristic length of various processes at the recombina- 
tion epoch reveals the physical process in the early uni- 
verse [2|. Length scales, which exhibit in the angular 
spectrum of the CMB as a characteristic feature, de- 
pend on the horizon length which is controlled by gravity. 
Usually, the horizon length at the recombination epoch 
is a function of the amount of non-relativistic matter. 
The strength of gravity caused by matter and radiation 
determines how much time consumed to make the uni- 
verse cool enough to recombine hydrogens. However, the 
strength of gravity depends not only on the amount of 
matter but also on the coupling strength of gravitational 
interaction. Thus the CMB anisotropy spectrum con- 
tains the information about the magnitude of the gravi- 
tational constant at the recombination epoch. 

Constancy of physical "constants" is the fundamen- 
tal issue which has long history and has attracted much 
interest Q. The recent attempts toward unifying all el- 
ementary forces Q predict the existence of scalar fields 
whose vacuum expectation values determine the physical 
"constants". In such context, scalar-tensor gravity the- 
ories, whose original version was proposed by Jordan Q 
and Brans and Dicke [sand were extended in a more 
general framework later |7j, are the most promising can- 
didates among the alternatives of Einstein gravity. In 
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scalar-tensor theories, the coupling of a massless scalar 
field to the Ricci scalar provides a natural framework of 
realizing the time- variation of the gravitational constant 
via the dynamics of the scalar field. 

In the Jordan-Brans-Dicke theory (hereafter, we 

refer to it as the Brans-Dicke theory for simplicity) which 
is the simplest example of scalar-tensor theories, a con- 
stant coupling parameter u> is introduced. In the limit 
lu — ► oo , the gravitational constant can not change and 
Einstein gravity is recovered. Although scalar-tensor the- 
ories including the Brans-Dicke theory are compatible 
with the Einstein gravity in several aspects, they have 
many deviations from it. Weak- field experimental tests 
in solar-system have constrained the post-Newtonian de- 
viation from the Einstein gravity, uj > 500 d, 0, 0, ^| . 
According to some recent reports, this bound would be 
updated to several thousands Such a small 

deviation implies that the variability of the background 
scalar field is also small. 

In more general scalar-tensor theories with non-trivial 
coupling functions w((f>), the small deviation is not the 
outcome of fine-tuning because the cosmological evolu- 
tion drives u> toward infinity in the late cosmological 
epochs and naturally reduces the present observable ef- 
fects of <f> field 15j. If such attractor mechanism takes 
place, the nature of gravity (the variability of the back- 
ground scalar field, weak-field deviations, etc.) can be 
significantly different in the early universe. Hence, in- 
formation on the different cosmological epochs may con- 
strain such theories. A simple and natural extension of 
the Brans-Dicke theory to the attractor model is the har- 
monic attractor model in which the scalar field has a 
quadratic effective potential of a positive curvature in 
the Einstein conformal frame. The analysis of big-bang 
nucleosynthesis (BBN) in the harmonic attractor model 
restricts two parameters characterizing the potential 
(its curvature /3 and today's gradient cyq). It is concluded 
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that the BBN limit on the possible deviation from Ein- 
stein gravity (2u>q + 3 = a ~ 2 ) is much stronger than the 
present observational limits in large j3(> 0.3) models. 
We can extract the information about the e arly universe 
also by the analysis of structure formation |l7j . The ad- 
vantage of the use of the CMB fluctuations is that the 
physics of CMB is well understood and that we now have 
very accurate observational data. The trace of primordial 
fluctuation can be seen clearly in the CMB anisotropy 
spectrum where the information on the early universe up 
to the last scattering time is projected on the acoustic 
peaks [3 |H H3- 

We compare the CMB temperature anisotropy spec- 
trum measured by WMAP with that in the scalar-tensor 
cosmological model mentioned above. The formulations 
of background and perturbation equations and the de- 
tailed explanation of physical processes can be found in 
our previous paper 20] . The purpose of this work is to 
constrain the deviation of the scalar-tensor gravity from 
the Einstein gravity and to find, in the framework of this 
model, how large the variation of the gravitational con- 
stant is allowed. 

This paper is organized as follows: In Section II, the 
scalar-tensor cosmological model and its observational 
consequences are explained. In Section III, we describe 
the overview of our analysis for constraining the scalar- 
tensor coupling parameters. In Section IV, the models 
are compared with the WMAP data and the result of \ 2 
test is presented. Finally, some conclusions are in Section 
V. 



where <f>o , &o and P are the present value of background 
4> field, today's potential gradient and curvature, respec- 
tively. We consider non-negative value of f3 since curva- 
ture near a minimum is positive. If f3 = 0, this model is 
reduced to the Brans-Dicke theory. Moreover, the model 
with «o — > and [3 = is the Einstein gravity. In the first 
post-Newtonian approximation, deviations from general 
relativity are proportional to the well-known Eddington 
(PPN) parameters as, 



7Edd - 1 = -2a 2 /(l + "o 2 ) , 
/?Edd-l = i/W/(l + ao 2 ) 2 . 



(3) 
(4) 



We see explicitly from Eqs. J3J and (J3J that post- 
Newtonian deviations from general relativity tend to zero 
with ao at least as fast as «q. This holds true for weak- 
field deviations of arbitrary post- Newtonian order [2*lj . 
One of the most stringent empirical limits for PPN pa- 
rameters is 

- 1.7 x 10~ 3 < 4 /? E dd - 7Edd - 3 < 1.5 x 10~ 4 (la) ,(5) 

which is obtained by the Lunar Laser Ranging experi- 
ment 22]. In the framework of the present model, this 
translates into 



a Q 2 < 1.5 x 10-7(/3 + l). 



(6) 



The cosmological evolution equations based on the the- 
ory are 



II. MODEL AND ITS PREDICTIONS 

The action describing a general massless scalar-tensor 
theory is 



S = 



1 



167rGn 



4>R 



(V0) 2 
+S m [ip,g^u], (1) 



where Go is the Newtonian gravitational constant mea- 
sured today, and w(</>) is the dimensionless coupling pa- 
rameter of the nonminimal coupling which is a function 
of (p. The last term in Eq.Q denotes the action of mat- 
ter which is a functional of the matter variable ip and the 
metric g^. The deviation from the Einstein gravity de- 
pends on the asymptotic value of <f> field at spatial infinity. 
According to the cosmological attractor scenario [13, the 
dynamics of the cosmological background 4> field is anal- 
ogous to that of a particle damping its motion toward the 
minimum of an external potential in the Einstein confor- 
mal frame. As the generic feature of a potential near a 
minimum is parabolic, we shall study the case where the 
potential is quadratic. This setup corresponds to w(0) of 
the following form, 



2w(0) + 3 = {a 2 - /?ln(0/<fo)}' 



(2) 
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{^G o a 2 (p-3 P )-0' 2 ^}. (9) 



A flat universe is assumed here. The prime denotes a 
derivative with respect to the conformal time, p and p are 
the total energy density and pressure, respectively. The 
effective gravitational constant measured by Cavendish- 
type experiments is given by [lOl l23| 



Gp 2cj(0)+4 
<t> 2u((j>) +3' 



(10) 



The requirement that today's gravitational constant be 
in agreement with the Newton's constant determines the 
present value of </> as 



4 + 2cj 
3 + 2w 



1 



ao 



(11) 

The system 



where c^o denotes the present value of uj(d} 
of perturbation equations is found in rcf. 20] . 

FigQ] shows the examples of typical <j> evolution. <f> is 
frozen during the radiation-dominated epoch and begins 
to grow at the matter-radiation equality time to real- 
ize the Newtonian gravitational constant at present. As 
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FIG. 1: The typical time evolution of <j> in the scalar-tensor 
models with ACDM parameters. The smaller <f> in the early 
epochs results in the smaller horizon length than that in the 
Einstein model. 



(f> moves upward, 2lu + 3 increases toward infinity, and 
therefore the present small deviation from the Einstein 
gravity is naturally realized. The increase of 2ui + 3 de- 
celerates the motion of cj) and finally <j) converges to some 
value which corresponds to the minimum of the effective 
potential. As increasing ao or /3, we obtain smaller ini- 
tial <j). This results in the smaller horizon length in the 
early epochs, which leaves distinct traces on the CMB 
spectrum. 

The typical CMB anisotropy spectra are shown in 
FigEl The locations of acoustic peaks (sound horizon 
scale) and that of diffusion damping tail are dependent 
on the horizon length at the recombination epoch. Since 
the matching condition of <po restricts the deviation of the 
present horizon length from that in the Einstein gravity, 
these angular scales directly represent the horizon length 
at recombination. Therefore, the shift of these angular 
scales to smaller scales is due to the large gravitational 
constant at recombination. The locations of acoustic 
peaks are shifted to smaller angular scales in proportion 
to the horizon length. Although the diffusion tail is also 
moved to smaller scale, it has weaker dependence on the 
horizon length, and hence the width between the first 
peak and the diffusion tail becomes thinner, which re- 
sults in the suppression of small scale peaks. The motion 
of <j> field and its fluctuation also distorts the spectrum. 
Roughly speaking, they make the first acoustic peak more 
prominent than higher peaks poj. However this effect is 
not so significant compared with the effect of changing 

"m0 or £!;,o- 



III. METHOD 

Let us describe the overview of our analysis. To con- 
strain the scalar-tensor coupling parameters, we compare 
the models to the WMAP temperature anisotropy spec- 



FIG. 2: The typical CMB temperature anisotropy spectra in 
the scalar-tensor models with ACDM parameters. Also the 
WMAP data is displayed for reference. The normalizations 
are marginalized to the WMAP data. The acoustic peak lo- 
cations are shifted to smaller angular scales due to the smaller 
sound horizon length (cf. Fig^l. 



trum 24] adopting the routine provided by the authors 
of ref . 25] to compute the values of % 2 . We employ a cos- 
mological constant (A) which is the simplest dark energy 
model that can account for the late time cosmic acceler- 
ation. 

In the model concerned in this paper, the CMB spec- 
trum depends on the three classes of model parame- 
ters. The first class consists of the conventional back- 
ground cosmological parameters which are today's Hub- 
ble parameter (/i), today's non-relativistic matter con- 
tent (Q m o), today's baryon content (Qbo)i the rcdshift 
of cosmic reionization (z re i ), today's CMB temperature 
(Tcmb), the helium mass fraction (Yff e ), and the neu- 
trino effective number (N„). We restrict ourselves to flat 
space models and hence today's amount of dark energy, 
which is introduced as a component of total energy and 
whose equation of state is pa/ PA = — 1> is a function of 
!] m o so that the modified Friedmann equation (Eq.©) 
is satisfied. The second class consists of perturbation pa- 
rameters characterizing the origin of fluctuations. They 
are the combination of adiabatic and isocurvature ini- 
tial perturbations, the scalar spectral index (n s ) and the 
overall normalization (A). The last class consists of the 
scalar-tensor coupling parameters which are the present 
deviation from the Einstein gravity (ao 2 ) and the curva- 
ture (/?) of the harmonic effective potential in the Ein- 
stein conformal frame. 

For the cosmological parameters, we use the follow- 
ing priors that are principally based on the WMAP 68% 
mean confidence range. Hubble parameter (h) takes 
the value between 0.67 and 0.77. As described in the 
next section, larger Q m o (smaller cosmological constant) 
shifts the acoustic peak locations to larger angular scales 
also in the scalar-tensor models. Therefore we allow 
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the non-relativistic matter of somewhat large content, 
!l m o £ (0.10,0.78), to confirm how small the value of 
X 2 for a set of (ao 2 ,(3) can be at its minimum. Sim- 
ilarly, the range of today's baryon content is set as 
ilbo £ (0.023,0.055). We employ instantaneous reioniza- 
tion at the redshift (z re io) between 12 and 22. We set the 
CMB temperature as T C mb = 2.72GK from COBE [U 
and the neutrino effective number (N v ) is fixed at 3.04. 
Although the helium mass fraction depends on the result 
of BBN, which is dependent on other parameters, it does 
not significantly affect the CMB spectrum. Therefore we 
set it as Y He = 0.24. 

According to the studies devoted to the generation 
of fluctuations durin g pa rticular inflation models with 
scalar-tensor gravity [27], |2^, |2(|, it is found that, al- 
though isocurvature perturbations could be produced 
during scalar-tensor inflation, they are in general neg- 
ligible compared with adiabatic perturbations and then 
the spectrum of the initial perturbations is not precisely 
scale invariant. Hence we employ adiabatic initial con- 
dition and allow the scalar spectral index to take the 
slightly deviated value from scale invariant one as fol- 
lows: n s £ (0.95, 1.03). As the overall normalization for 
each model, we survey the region between 0.8 and 2.2 in 
unit of COBE normalization factor. 

Since we are interested in setting a constraint from 
CMB alone, we survey the ranges of ao and (3 includ- 
ing the regions which have been ruled out by the solar 
constraint shown in Eq.JBJ as 

a 2 E (4 x 10~ 8 ,4 x 10~ 2 ), (12) 
/3 £ (0,0.45). (13) 

These regions roughly correspond to 10 < u>q < 10 from 
Eq.@. 



IV. COMPARISON WITH OBSERVATIONS 

In this section, we show the result of \ 2 test comparing 
the harmonic attractor model with the WMAP data. 

In FigOU we show the \ 2 contour map on ao 2 — (3 plane, 
marginalizing over the other parameters. We find that 
the scalar-tensor coupling parameters are constrained as 

a 2 < 5 x 10- 4 - 7/3 (2a), (14) 
a 2 < 10- 2 - 713 (4(i). (15) 

The contour map in Fig|3|has a sharp edge approximately 
on the curve, a 2 = 10~ 2 ~ 7/3 . Beyond this curve, x 2 of the 
models rapidly increases and then those models are sta- 
tistically improbable to explain the observed spectrum. 

Next, we consider the variation of the gravitational 
constant at the recombination epoch. In Fig0J we show 
the x 2 for the gravitational constant at the recombination 
epoch (G rec = G(4> rec )), marginalizing over other param- 
eters. Here, <p rec is the value of <j) at the recombination 
epoch. We find that the deviation of the gravitational 



constant from today's value is constrained as 

|G rec -G |/G <0.05 (2a), (16) 
|G rec -G |/G <0.23 (4a). (17) 

The deviation of the gravitational constant from today's 
value up to 5 % does not significantly change \ 2 an d 
the degree of fit of such models is comparable to that of 
WMAP team's best fit model, while even larger deviation 
increases % 2 . In the models on the curve, the locations 
of acoustic peaks are fitted to those of the observed spec- 
trum and then, in order to pull back the peaks shifted by 
the larger gravitational constant to observed locations, 
f2 m o and fi&o are out of their favorable values to fit the 
shape of the spectrum. The x 2 for the gravitational con- 
stant at our initial time (Gini = G(<j>i n ij) is also displayed 
for reference. Here, 4>ini is the value of 4> at Zj n j(~ 10 8 ). 
We find that the CMB constraint on the gravitational 
constant at z, n j, \Gj n j — Gpl/Go < 0.12, is on the same or- 
der as the BBN bound 30]: 0.7 < G B bn/G q < 1.4(2a), 
and it can be comparable to the BBN bound on the har- 
monic model Eq.J2l 16j since our analysis is limited to 
13 < 0.45. 

In Fig|Sl we show the allowed post-Newtonian devia- 
tions at three different epochs. The curve for today's de- 
viation parameter (ao 2 ) is identical to the cross section of 
the contour map at (3 = and hence it is identical to the 
X 2 curve for the Brans-Dicke models. The sharp edge 
mentioned above is located around ag 2 ~ 10 -2 (a->o ~ 
50). The models whose Brans-Dicke parameter is much 
smaller than the lower bound by the solar system exper- 
iments can be compatible with CMB fluctuation, which 
was expected previously. Even if we require that the de- 
gree of fit should be comparable to that of WMAP team's 
best fit model, the boundary of the corresponding region 
is still beyond the solar bound. On the other hand, the 
deviation at the recombination epoch and Zi n i are rela- 
tively loosely constrained: a rec 2 — l/(2u(<f) rec ) + 3) < 
7 x 10- 2 , a ml 2 = l/(2cj(0 ini ) + 3) < 2 x lO" 1 ^). This 
mig ht have implications for extended inflation scenarios 

In FigsEl and [7| we plot the \ 2 f° r ^mo and flbo, 
marginalizing over other parameters. Compared with 
the Einstein models, the scalar-tensor models can be 
more probable in high density models because the scalar- 
tensor models have the two more tunable parameters to 
fit acoustic peak locations. Although the peak locations 
are dependent on baryon abundance, large amount of 
non-relativistic matter results in very low peak heights 
which cannot be compensated by baryon drag especially 
on the 2nd peak. On the other hand, in low £lbo models, 
there is no definite difference between the scalar-tensor 
and the Einstein models. The observed spectrum have 
the prominent peaks and hence low density and small 
baryon abundance models arc not so ill-fitted compared 
with the scalar-tensor models within the surveyed param- 
eter range. 

Finally, we comment on the case if we allow non-flat 
universe models. Even in the case of non-flat models, 
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FIG. 3: The x 2 contour map on ao 2 — f3 plane for the 
scalar-tensor ACDM models, where the other parameters are 
marginalized. Also the solar bound from the LLR experiment 
is shown for reference. Note that the horizontal axis ao 2 is 
almost proportional to yEdd — 1 from Eq.©. There exists a 
sharp edge approximately on the curve, a 2 , = 10~ 2 ~ 713 . Be- 
yond the curve, x 2 of the models rapidly increases. Although, 
in the figure, the boundary curves only up to 4<r level are 
drawn, % 2 continues to inflate for larger coupling parameters. 

the constraint for the scalar-tensor coupling parameters 
(ao 2 ,P) of Ea. (fH))l would not significantly change be- 
cause the value of <j> in the early epochs depends on these 
parameters exponentially. On the other hand, the con- 
straint for the gravitational constant at the recombina- 
tion epoch would become much weaker because acoustic 
peak locations can be easily modulated in curved models 
without disturbing CDM or baryon abundance signifi- 
cantly. This degeneracy would be improved by obser- 
vations at even smaller scales uncovering the diffusion 
cut-off IH because the diffusion length depends on hori- 
zon length in a different way from peak locations. The 
ratio of the angular scale of the sound horizon to that 
of the diffusion scale is independent of projection effect, 
and hence it can provide the information on the horizon 
length at the recombination epoch. 



V. CONCLUSION 

In this paper, we have quantitatively compared the 
CMB temperature spectrum in the scalar-tensor flat 
ACDM model to the WMAP data by x 2 test of good- 
ness of fit. 

The present deviation from the Einstein gravity (ckq) 
must be smaller than 5 x l(T 4 - 7/3 (2a), and 10~ 2 ~ 7/3 
(4<r) for < < 0.45. This constraint is much stronger 
than the bound from the solar system experiments for 
large (3 models, i.e., (3 > 0.2 and 0.3 in 2a and 4cr lim- 
its, respectively. Within the framework of the harmonic 
attractor model, the difference between the gravitational 
constant at the recombination epoch and at the present 
is constrained as \G(z = z rec ) — Gq\/Gq < 0.05(2<r) and 
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FIG. 4: The \ 2 f° r the scalar-tensor ACDM models as a 
function of G rec and d n i, where the other parameters are 
marginalized. The degree of freedom v{= cr 2 /2) is 891. 
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FIG. 5: The x 2 f° r the scalar-tensor ACDM models as a 
function of ao 2 ,a rec 2 and cti„i 2 , where the other parameters 
are marginalized. The curve for ao 2 is idential to the % 2 f° r 
the Brans-Dicke models. The degree of freedom v{= cr 2 /2) is 
891. 

0.23(4cr). This is the first-time bound on the variation of 
the gravitational constant from CMB anisotropy spec- 
trum. While the present deviation from the Einstein 
gravity is severely constrained (a 2 , < 10~ 4 ), larger de- 
viation during radiation dominated epochs is compatible 
with CMB. Indeed a™ 2 up to 2 x 10~ 2 can be within 2a 
level. Although our analysis is limited to flat models, the 
further small scale CMB anisotropy data which will be 
provided in future |35| would break the degeneracy and 
would significantly improve the bound. 
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